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ABSTRACT 

We present high angular resolution VLA and BIMA observations of NH3, 
HGO+, HCN and SO molecular emission and 1.4, 3.5 and 7 mm continuum 
emission toward the exciting source of the HH 80-81 system. This object is one 
of the few massive protostars known to be driving a coUimated outflow. We 
report the first detection of SO 85-44 molecular emission toward the exciting 
source of HH 80-81, suggesting that this transition may be a good tracer of 
molecular gas near massive protostars. We also detected toward this source dust 
continuum emission at 1.4 and 3.5 mm. From the SO molecular emission and 
the dust emission we roughly estimated that the molecular mass associated with 
the circumstellar surroundings of the exciting source of the thermal jet is in the 
range 1-3 M©. Weak and broad (2,2) ammonia emission was also found in the 
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direction of the jet suggesting the presence of small amounts of molecular gas at 
high temperatures (> 50 K). 

The VLA observations show the presence of three ammonia components to- 
ward the HH 80-81 region. The brightest component peaks at ~8" to the NE 
of the thermal jet and is associated with the H2O maser spots in the region. A 
second ammonia clump is located about 25" to the NE of the jet and is asso- 
ciated with class I methanol masers. The third ammonia component is located 
1' to the south from the thermal jet and may be a molecular core yet without 
stellar formation. The BIMA observations show that the strongest emission in 
the HCO''" and HCN hnes originate close to the H2O maser, and cover the same 
spatial region and velocity range as the brightest ammonia component. 

Subject headings: ISM: outflows— individual (HH 80-81, GGD 27, IRAS 18162- 
2048) — ISM: HH regions — molecules — radio lines — ISM: Herbig-Haro objects — 
astrochemistry 

1. Introduction 

The HH 80-81 system, also known as GGD 27 (Gyulbudaghian ct al. 1978), consists 
of a highly collimatcd radio jet that powers an outflow with a total projected extent of 5.3 
pc, at an adopted distance of 1.7 kpc (Rodriguez et al. 1980; Martf, Rodriguez, & Reipurth 
1993; 1995; 1998). The thermal radio jet associated with HH 80-81 terminates to the south 
in the extremely bright optical HH objects 80 and 81 (Reipurth & Graham 1988; Marti, 
Rodriguez, & Reipurth 1993, 1995; Heathcote, Reipurth, & Raga 1998) and to the north in 
the radio source HH 80 North (Marti et al. 1993; Girart et al. 1994; 2001). The central radio 
continuum source has a bright far-infrared counterpart, IRAS 18162—2048, which implies 
the presence of a luminous (~ 2 xlO"^ L0) young star or cluster of stars (Aspin & Geballe 
1992). Water maser emission has been reported in the vicinity of the central radio continuum 
source (Rodriguez et al. 1978; 1980). The H2O masers do not coincide with the jet core, 
suggesting that they are being powered by a different star (Gomez, Rodriguez, & Marti 1995; 
Martf, Rodriguez, & Torrelles 1999). Methanol emission at 95 GHz has been reported by 
Val'tts et al. (2000) to the northeast of the exciting source of the jet at the position of 
q;(2000)= 18'^ 19"* 12f89; 5(2000)= -20° 47' 13'.'6 with a positional error of 10". These 
methanol masers are classified as class I which are believed to trace high-mass protostars 
that have not significantly ionized their environment. 

The most remarkable characteristic of the HH 80-81 system is the high degree of col- 
limation of the jet, with an opening angle of ~1° (Martf, Rodriguez, & Reipurth 1993). 
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CoUimated outflows and accretion disks are thought to be present at the same evolutionary 
stage in the formation of both low-mass and high-mass stars (e.g. Garay & Lizano 1999; 
Kurtz 2000). An accretion disk around the central star is believed to be responsible for the 
ejection and collimation of the gas. In contrast with low-mass stars, where disk structures 
have been observed commonly around young objects (Beckwith & Sargent 1996; Wilner, Ho, 
& Rodriguez 1996; Rodriguez, et al. 1998; Wilner et al. 2000), in high- mass stars it has 
been difficult to establish the presence of accretion disks. Whether the massive protostars 
are formed by accretion (as it is the case for low mass stars) or by merging of protostellar 
condensations (Bonnell & Bate 2002) remains controversial. It has been argued, for example, 
that the ionizing flux present in massive stars could rapidly disrupt the disk. However, the 
presence of well coUimated outflows in a few objects also points to the presence of accretion 
disks surrounding massive stars. If massive OB stars are formed by accretion, we expect 
that disks and jets will be present in their earliest stages of evolution, as in the case of low 
mass stars. Accretion disks and collimated outflows around B-type protostars have been 
tentatively reported toward Cepheus A HW2 (Rodriguez et al. 1994; Torrellcs et al. 1996), 
IRAS 20126+4104 (Cesaroni et al. 1997; Zhang, Hunter & Sridharan 1998), G192.16-3.82 
(Shepherd et al. 1998; Devine et al. 1999; Shepherd, Claussen, & Kurtz 2001), AFGL 5142 
(Zhang et al. 2002) and AFGL 2591 (Trinidad et al. 2003). A possible coUimated outflow 
associated with a luminous 0-type protostar is IRAS 16547—4247, a bright infrared source 
that is associated with a triple radio source consisting of a compact central object and two 
outer lobes (Garay et al. 2003). There is, however, no known unquestionable case of a 
circumstellar disk associated with a massive 0-type protostar. 

The luminosity of HH 80-81 is ~2xl0^ L©, equivalent to that of a BO ZAMS star. Since 
HH 80-81 possesses an extraordinary well coUimated jet, it seems reasonable to propose that 
it should be originating from an accretion disk. However, until now there has been no clear 
evidence for the existence of this disk. A bipolar CO outflow has been detected and mapped 
in the direction of HH 80-81 by Yamashita et al. (1989) which spatially coincides with the 
infrared reflection nebula in the region. Dense molecular gas traced by NH3 (1,1) emission 
was reported by Torrelles et al. (1986) from single dish observations, whereas Girart et al. 
(1994) reported the detection of (1,1) and (2,2) inversion transitions of ammonia downstream 
in HH 80 North. Also a CS molecular clump has been observed towards the exciting source 
of HH 80-81 which exhibits an elongated structure (~60 x 25") perpendicular to the CO 
molecular outflow, that has been interpreted as a molecular disk (Yamashita et al. 1989). 
However, high angular resolution observations are needed in order to establish the existence 
of a more compact circumstellar disk. We consider HH 80-81 as one of the best candidates 
for the detection of a protoplanetary disk using molecular lines. Accordingly, the main goal 
of this work was to make observations, in several molecular lines and with adequate angular 
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resolution, towards the core of HH 80-81 to search for an accretion disk that could account 
for the coUimated outflow. 

2. Observations 
2.1. VLA 

The 1.3 cm line and continuum observations were made on 1999 April 12 with the Very 
Large Array (VLA) of the NRAO^ . The array was in the D configuration giving an angular 
resolution of 4" at 1.3 cm. Structures larger than 60" were undetectable. We observed the 
(J,K) = (1,1) and (2,2) inversion transitions of the ammonia molecule in the spectral line 
mode towards the core of the HH 80-81 complex. We assumed rest frequencies of 23694.496 
and 23722.633 MHz for the (1,1) and (2,2) transitions, respectively. We used a bandpass 
of 3.125 MHz, centered at an LSR velocity of 10 km s~^ and 63 spectral channels 48.8 kHz 
wide each (~ 0.6 km s~^ at the observing frequencies) plus a continuum channel containing 
the central 75% of the total band. Flux densities were bootstrapped from observations of 
1331-1-305, which was assumed to have a flux density of 2.4 Jy at 1.3 cm. The bandpass 
calibrator was 1229-1-020 with a flux density of 28.6 Jy and the phase calibrator was 1733-130 
with a flux density of 3.8 Jy at 1.3 cm. The data were edited and calibrated in the standard 
manner using the Astronomical Image Processing System (AlPS) developed by the NRAO. 

All spectra were Hanning smoothed to a velocity resolution of 1.2 km s^^. The rms noise 
levels in a single spectral channel, after Hanning smoothing, were 2.5 and 2.2 mJy beam~^ 
(corresponding to 0.38 and 0.33 K) for the (1,1) and (2,2) observations, respectively. The 
system temperature during observations was in the order of 150 K. We applied a primary 
beam correction (with the task PBCOR of AIPS) to the ammonia data in order to have a 
good estimate of the flux for the emission that is displaced from the phase center. 



^The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooper- 
ative agreement with the National Science Foundation 
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2.2. 



BIMA 



The observations were carried out in 1999 September 27 at 3.5 mm and 2001 May 6 
at 1.4 mm with the 10-antenna BIMA array^ in the C configuration. The phase cahbrator 
used was 1733-130. Absolute fiux calibration was performed by observing Mars. System 
temperatures were in the 200-400 K range at 3.5 mm and in the 1600-4000 K range at 
1.4 mm. The 3.5 mm frequency setup included the HCO"*" 1-0 and HCN 1-0 lines, whereas 
the 1.4 mm setup included the SO 55-44 and H2CO 3o,3-2o,2 hues. Self-calibration was 
performed in the 1.4 mm data by using the continuum emission. Maps were made with the 
{u, v) data weighted by the associated system temperatures. More details of the observations 
are given in Table 1. 



Using the VLA we detected an unresolved (< 2''6) 1.3 cm continuum source coincident 
in position with the exciting source of the thermal jet. The estimated fiux density at 1.3 cm 
is 6 ± 1 mJy. This is consistent with the value reported by Marti, Rodriguez & Torrelles 
(1999). The 7 mm flux density was obtained from data taken with the VLA on 1994 April 
15. The unresolved (< 0'.'4) 7 mm source is shown in Figure 1. The continuum images at 3.5 
and 1.4 mm made with BIMA are shown on Figure 2. To discuss the nature of this emission 
we present in Figure 3 the continuum spectrum towards the exciting source of the thermal 
jet at cm and mm wavelengths. The 20, 6, and 2 cm fiux densities are taken from Marti et 
al. (1993), while the 1.3 cm data point is that presented above from our VLA observations. 
The fiux densities at 7, 3.5 and 1.4 mm toward the exciting source of the outfiow are listed 
in Table 2. From Figure 3 we can see that the cm emission shows the characteristic fiat 
spectrum of nearly opticahy thin free-free emission (Si, ~ i/'^-^='=°-^), but that at 7 mm and 
shorter wavelengths, the flux density rises rapidly. This spectrum suggests that the flux 
densities at 3.5 and 1.4 mm are completely dominated by the dust (at 3.5 mm the expected 
contribution from the thermal jet is only ~0.01 Jy). The spectral index derived from the 
7, 3.5 and 1.4 mm fluxes is 2.6, consistent with optically-thin dust emission. There is also 
1.4 mm emission associated with the C ammonia component (see Figure 2, and Table 3). 



^Thc BIMA array is operated by the Berkeley- Illinois-Maryland Association with support from the Na- 
tional Science Foundation. 
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RESULTS 



3.1. 



VLA and BIMA continuum emission 
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3.2. VLA spectral line emission 

The (1,1) and (2,2) main inversion transitions of ammonia were strongly detected in 
emission toward the northeast of the jet, but only very weak (2,2) emission was observed at 
the position of the exciting source of the thermal jet. 

Figure 4 shows the velocity-integrated ammonia (1,1) emission toward HH 80-81. We 
identify three main ammonia structures; the bright, C component which exhibits an extended 
structure and peaks close to the position of the water masers in the region, a small ammonia 
clump to the northeast of the C component (hereinafter called the NE component) and the 
S component which is located ~1' to the south of the C component. 

Figure 5 shows the channel maps of the NH3 (1,1) and (2,2) emissions toward the 
surroundings of the exciting source of HH 80-81. In general, the ammonia emission from the 
molecular cloud in both transitions appears distributed to the northeast of the thermal jet, 
spanning over the same velocity range, from ~ 10 to 14 km s~^. As noted before, the peak 
emission of the C ammonia component is nearly coincident with the position of the H2O 
maser spots (Gomez, Marti, & Rodriguez 1995; Marti, Rodriguez & Torrefies 1999). 

Figure 6 presents a close-up of the velocity-integrated emission in the (1,1) and (2,2) 
ammonia lines from the C component and the (1,1) line from the S component. The top 
panel also shows the (1,1) ammonia emission from the NE ammonia clump where methanol 
maser emission has been reported (Val'tts et al. 2000). The S component was detected 
only in the (1,1) ammonia transition and in the velocity range from 11.9 to 13.7 km s~^, 
suggesting it corresponds to a cold molecular clump. 

The spectra of the NH3 (1,1) and (2,2) inversion transitions from each individual am- 
monia component (integrated over their whole area of emission) are presented in Figure 7. 
The (1,1) main and satellite lines were detected toward each individual ammonia compo- 
nent, while in the (2,2) ammonia transition emission was detected only in the main line 
and toward the C and NE components. In Table 3 we list the ammonia parameters of the 
main and satellite ammonia lines from the C, NE and S components determined by fitting a 
Gaussian profile to the corresponding spectrum. The peak intensity of these three ammonia 
components are around the LSR velocity of -|-12 km s~^. 

3.3. BIMA spectral line emission 

Figure 8 shows line profiles and maps of the SO 85-44, HCO"'" 1-0 and HCN 1-0 emission. 
The SO emission arises from a compact region around the thermal jet. A Gaussian fit to 
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the SO profile gives a peak flux density of 7.6 ± 1.6 Jy, a central velocity of Ylsr— 12.1 
± 0.2 km and a linewidth of AV= 1.3 ± 0.4 km s~^. This is the first time that SO 
emission is detected toward the exciting source of HH 80-81 and suggests that this transition 
may be a good tracer of molecular gas near massive protostars. The HCO+ and HON 
emissions trace basically the same spatial region as the ammonia central component and in 
a similar velocity range. Weak HCO"'" emission is observed toward the NE ammonia clump, 
confirming its existence. There is some enhancement of the HCO"*" and HON emission toward 
the direction of the thermal radio jet, but it is difficult to separate the molecular gas arising 
from the jet region from the more extended central component. Despite the similitude of the 
HCO"*" and HON integrated maps with the integrated ammonia maps, the observed HCO+ 
and HCN spectra show a clear self-absorption feature at the center of their lines (see Fig. 8), 
which is not observed in the ammonia lines. It is possibly a true self-absorption since the 
shortest baselines also show this absorption feature. 

We also searched for emission in the H2CO 3o,3-2o,2 line toward the region. The inte- 
grated velocity spectrum in a velocity range from 10 to 14 km s~^ did not show emission at 
a 3-(T level of ~2 Jy beam~^. 

4. DISCUSSION 

The line data presented here reveal the presence of three main molecular structures that 
we called the C, the NE, and the S components. In what follows we discuss each of these 
components and derive their bulk physical parameters, and conclude with a discussion of 
the molecular and dust emission toward the surroundings of the jet. 

4.1. The C component 

This large structure of molecular gas has a peak position located about 8" northeast 
from the exciting source of the thermal jet. Emission toward this component was detected 
in the (1,1) and (2,2) transitions of NH3, and HCN (1-0) and HCO+ (1-0) lines, indicating 
the presence of considerable amounts of cold molecular gas. The strongest H2O maser in the 
region (marked with a triangle in Figures 4 to 8; with a Vlsr = ~61 km s^^ and a peak 
position at a (2000)= 18^* 19™ 12f509; 5(2000)= -20° 47' 27!'42, Gomez et al. 1995), lies close 
(2") to the molecular peak emission of the C component (see Table 3). The LSR velocities 
of the individual H2O maser spots are between —50 to —72 km s~^ (Marti et al. 1999), 
considerably blueshifted with respect to the NH3, HCN and HCO"*" molecular emission (-1-12 
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km s A deep 3.5 cm radio continuum map reveals the presence of a weak (~0.05 mJy) 
unresolved (< 0''2) radio source (VLA 3 in Gomez et al. 1995; a(2000)= 18'^ 19"^ 12f48; 
5(2000)= —20° 47' 27''42) that is very close to the strongest H2O maser. The position of 
VLA 3 is also located a few arcsec from the position of the 1.4 mm emission detected by us 
(a(2000)= 18^ 19™ 12f 13; 5(2000)= -20° 47' 30''57; see Fig. 2) and we tentatively propose 
that they correspond to the same source. Assuming that the radio continuum emission from 
VLA 3 arises from an HII region, it requires a B3 ZAMS star to maintain its ionization. All 
the above data support the idea that the C molecular clump harbors a cluster of stars in an 
early stage of formation. Mid-IR observations indicate the presence of a pre-main sequence 
stellar cluster (Aspin & Geballe 1992; Yamashita et al. 1987) at this position. 

The derived parameters of the C ammonia component, as well as of the other ammonia 
components, are summarized in Table 4. The C ammonia source is the warmer component 
having a rotational temperature of ~ 24 K. It shows a velocity gradient of ~3 km s~^ over 
an angular scale of ~12" (see Fig. 9). If due to rotation, it implies a mass of ~30 Mq, 
which is in good agreement with the mass of ~40 Mq estimated from the virial assumption 
(Table 4). 

4.2. The NE component 

The emission from the NE ammonia clump is seen in the velocity range from 10.0 to 
10.6 km s"-*^. Neither HCN nor radio continuum emission was detected toward this clump. 
This component is associated with methanol masers (coincident within the positional error 
of the masers of 10"; see Figure 3) as well as with weak HCO+ emission. The 95 GHz class I 
methanol masers found are in a velocity range from 12.4 to 13.7 km s~^ (Val'tts et al. 2000), 
slightly redshifted with respect to the LSR velocity of the ammonia clump (~11 km s~^). 
Class I methanol masers are characterized by a lack of radio continuum emission, and may 
be related with the first stages of massive stellar evolution (Garay & Lizano 1999). The 
relatively warm rotational temperature of this component (22 K; see Table 4) also favors 
the presence of embedded stars. The molecular mass estimated from the ammonia data 
is ~ 1 Mq. On the other hand, the virial mass is aproximately 13 Mq, suggesting that 
this molecular clump may not be in equilibrium. Most probably, this warm clump is a 
substructure of the C component. In fact, from the channel maps of the (1,1) ammonia 
emission (Fig. 5) we note that at 11.9 km s^^ there arc several weak extended structures 
around the C component. In this channel map we also observe an EW filament resembling 
the "ammonia fingers" in OMC-1 (Wiseman & Ho 1998). 



-9- 



4.3. The S component 

This molecular clump was detected only in the (1,1) main transition of NH3 and ap- 
pears located toward the south of the exciting source of the jet. We did not detect emission 
of HCO+ and HCO at a 3a level of 1.5 Jy beam~^. The ammonia (1,1) line has Vlsr ^ 
13 km s~^, slightly redshifted with respect to the C component. The low rotational tem- 
perature of this component (10 K; see Table 4) suggests that it does not have embedded an 
object of large luminosity. The estimates of the molecular mass and the virial mass are very 
similar, ~ 13 M0, suggesting that this clump is probably in equilibrium. It will then be 
interesting to study this component in detail since it may be a precursor to the formation of 
stars. Marti et al. (1993) found that along the north-south direction of the thermal jet there 
are several knots that emit at radio frequencies and in particular the S ammonia component 
is close to one of those knots (source VLA 13 from Marti et al. 1993): 21" or 0.18 pc in 
projection. No counterpart was found at other frequencies for this new S ammonia clump 
using the SIMBAD database. An inspection of the 2MASS All Sky Survey shows that there 
are many near-IR sources in the GGD 27 region. Yet, there is a region of roughly 40" of 
diameter around the S ammonia clump where there are no near-IR sources. This suggests 
that the S ammonia clump is part of a very dark molecular cloud, where apparently there 
is no star formation going on. The presence of a possible shock excited radio source nearby 
(source VLA 13) suggests that this clump may belong to the class of irradiated clumps, that 
is, clumps that are exposed to the strong radiation of a Herbig-Haro object (e.g. Girart et al. 
1994, 2002). In these clumps the ammonia emission is strongly enhanced (Viti & Williams 
1999; Viti et al. 2003). If the S clump is an irradiated clump, the estimated mass would be 
much lower than 13 Mq. 



4.4. Emission Toward the Jet 

Our results show the presence of SO molecular emission and continuum emission at 
mm wavelengths in the direction of the exciting source of the HH 80-81 jet. The 65-44 
transition of SO is a tracer of warm (>40 K) and dense (~10^ cm~^) gas (Chernin et al. 
1994). The SO emission is detected at > 4(j level in both the spectrum and the integrated 
map of Figure 8. Prom the SO line emission and assuming optically thin conditions, we 
obtain a beam-averaged column density of N(SO) in the range 2.5 — 3.4 x 10^^ cm~^ for a 
temperature range of 50 to 100 K. Assuming that the SO/H2 ratio is in the 0.5 — 2.0 x 10~^ 
range (Helmich & van Dishoeck 1997), we obtain a beam-averaged H2 column density in the 
range of 1.3 — 6.8 x 10^^ cm~^, that translates in a mass of 0.7 — 3.5 Mq. On the other hand, 
assuming that the dust emission is optically thin and that the dust temperature is in the 50 
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— 100 K range, we derive a total mass from the dust emission of 3.7 — 1.8 Mq (Beckwith 
& Sargent 1991). We conclude that both the SO and the dust data are consistent with the 
presence of a few solar masses (1 — 3 Mq) of gas closely associated with the surroundings 
of the exciting source of the jet. 

We found that no strong ammonia emission was detected toward the jet. The position 
of the exciting source of the thermal jet lies at the edge of the C ammonia component (see 
Figure 4). In the image of the C component (Figure 6) we see a small enhancement of the 
(2,2) ammonia emission along the direction of the thermal jet, while there is no emission 
in the (1,1) image. Figure 10 shows the spectra of the (1,1) and (2,2) inversion transitions 
at the peak position of the thermal jet. We detect a weak and broad (2,2) ammonia line 
emission which suggests that warm ammonia gas is present in the vicinity of the exciting 
source of the thermal jet. From a Gaussian profile fit to the (2,2) spectrum we obtain a peak 
flux of 4.8 ± 0.7 mJy beam~^, an LSR line center velocity of 11.5 ± 0.3 km s~^ and a line 
width of 5.4 ± 0.8 km s~^. Using these parameters for the (2,2) line and a 2a upper limit 
for the (1,1) ammonia line of ~3 mJy beam~^, it is possible to estimate a lower limit to the 
rotational temperature of ~50 K. This is unusually warm molecular gas, that we suggest is 
being heated by the star that produces the jet. However, the determination of its precise 
location and nature cannot be addressed with the data presented here. A crude upper limit 
for the molecular mass associated with the jet can be obtained as follows. From the peak 
(2,2) hue flux density of 4.8 mJy beam~^ and the angular resolution of the observations, we 
estimate a line brightness temperature of ~0.7 K. Assuming that the excitation temperature 
is 50 K,this implies r(2,2,m) ~ 0.014, an optically thin line. Adopting a hnc width of 5.4 
km s~^, we derive a total column density of ammonia of N^NH^) ~ 2 x 10^^ cm~^. Using the 
ratio of [/J2/NH3] = 10^ used for the ammonia components, this implies N{H2) ^ 2 x 10^^ 
cm^^. Finally, assuming that the gas is distributed over a physical size similar to the beam, 
we obtain a mass estimate of 0.3 Mq. It is unclear if this molecular gas is associated with 
the jet or with the edge of the C molecular component. In any case, the mass derived from 
ammonia for the surroundings of the jet is an order of magnitude smaller than that derived 
from SO or dust, suggesting that ammonia is strongly depleted. 



5. Conclusions 

We have made high angular resolution VLA and BIMA observations of NII3, HCO^, 
HCN, and SO molecular gas as well as 1.4, 3.5 and 7 mm continuum emission toward the 
exciting source of HH 80-81. Our main conclusions are summarized as follows. 

We detected millimeter continuum emission (1.4, 3.5 and 7 mm), as well as SO 55-44 



- 11 - 



line emission, associated with the exciting source of the HH 80-81 outflow. These emissions 
may be arising from a circumstellar disk, but higher angular resolution data are required to 
study the morphology of the emitting dust and gas. This is the first time that SO 65-44 line 
emission is detected toward the exciting source of HH 80-81, suggesting that this transition 
may be a good tracer of molecular gas near massive protostars. Only weak ammonia (2,2) 
emission was detected toward the exciting source of the jet, indicating a temperature ~50 
K. Unfortunately, the precise location and nature of the warm molecular gas traced by the 
ammonia cannot be determined with the present data. 

Three molecular components have been detected in the surroundings of the exciting 
source of the HH 80-81 outflow. The brightest of these condensations (the C component) was 
detected in NH3, HCO"*", HCN and appears associated with H2O maser emission suggesting 
that this component harbors a cluster of young stars. The NE component was detected 
in NH3 and HCO"^ and is associated with class I methanol masers. The S component was 
detected only in NH3, and seems to be a massive (about 15 M0) and cold clump without on- 
going stellar formation. Alternatively, this component may belong to the class of irradiated 
clumps, with VLA 13 being the source of the radiation. 

YG and LFR acknowledge financial support from DGAPA-UNAM and CONACyT, 
Mexico. JMG is supported by MCyT grant AYA2002-00205. GG acknowledges support 
from FONDECYT project 1010531. JM also acknowledges partial support by DGI of the 
Ministerio de Ciencia y Tecnologia (Spain) under grant AYA2001-3092, as well as by Junta 
dc Andalucia (Spain) under PAl group FQM322. This research has made use of the SIM- 
BAD database, operated at CDS, Strasbourg, France. The 2MASS project is a collabora- 
tion between The University of Massachusetts and the infrared processing Analysis Center 
(JPL/Caltech), with funding provided primarily by NASA and NSF. 
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Table 1. PARAMETERS OF THE BIMA OBSERVATIONS 







1/ 


Synthesized Beam 


Bandwidth 




RMS^ 


Epoch 


Observation 


(GHz) 


HPBW 


RA. 


(MHz) 


(kms-i) 


(Jy Beam"-^) 


1999 Sep 27 


Continuum 


85.12 


19'.'0 X 6'.'3 


-14° 


100 




0.021 




HCO+ 1-0 


89.18852 


18'.'7 X 6'.'1 


-14° 


25 


0.66 


0.56 




HON 1-0 


88.63185 


18'.'3 X 6'/l 


-14° 


25 


0.66 


0.52 


2001 May 6 


Continuum'^ 


216.77 


TA X 2'/9 


-7° 


1200 




0.038 




SO 55-44 


215.22065 


TA X 3'/l 


-9° 


25 


0.68 


1.80 




H2CO 3o,3-2o,2 


218.22219 


TA X 3'/l 


-9° 


25 


0.67 


2.00 



^Spectral resolution. 

^Pcr channel. The brightness temperature of the 3.4 mm data is Tb[K) ~ 1.4 x 10~^5'(mJy), while 
for the 1.4 mm data is Tb{K) ^ 1.1 x lQ-^S{mJy). 

^The continuum maps at 1.4 mm were made using the whole upper and lower sidebands. 



Table 2. DUST EMISSION 



A Total Flux (Jy) 

(mm) Thermal Jet C com])oiient 



7.0 0.014±0.002 < 0.006 

3.5 0.10±0.02 < 0.08 

1.4 0.97±0.11 0.21±0.05 
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Table 3. AMMONIA (1,1) AND (2,2) LINE PARAMETERS 



Peak Position 

Component a(2000) 5(2000) Vlsr AV Flux Density" 6^ 

(h m s) (o / //) (km s"^) (km s"^) (mJy) (//) 



c 


(l,l;m) 


18 


19 


12.54 


-20 47 25.5 


12.2 


± 


0.1 


2.3 ± 0.1 


260 ± 6 


8.9 




(i,i;s) 










12.1 


± 


0.1 


2.6 ± 0.1 


80 ± 6 






(Ll;s) 










12.2 


± 


0.1 


1.9 ± 0.1 


108 ± 6 




c 


(2,2;m) 


18 


19 


12.55 


-20 47 25.4 


12.2 




0.1 


2.4 ± 0.1 


188 ± 6 


10.6 


NE 


(l,l;m) 


18 


19 


13.52 


-20 47 13.2 


10.6 


± 


0.1 


2.1 ± 0.2 


25 ± 2 


3.5 




(i,i;s) 










11.7 


± 


0.4 


3.4 ±0.8 


8 ± 2 






(i,i;s) 










11.3 


± 


0.4 


3.9 ±0.9 


8 ± 2 




NE 


(2,2;m) 


18 


19 


13.42 


-20 47 13.1 


10.5 


± 


0.1 


1.9 ± 0.3 


14 ± 2 




S 


(l,l;m) 


18 


19 


11.94 


-20 48 29.3 


13.0 


± 


0.1 


1.6 ± 0.1 


111 ± 6 


6.1 




(i,i;s) 










12.9 


± 


0.1 


1.6 ± 0.2 


67 ± 6 






(i,i;s) 










13.0 


± 


0.1 


1.7 ± 0.2 


58 ± 6 




S 


(2,2;m) 


















< 18 





^At the peak of the line. 

''Angular size estimated from a Gaussian fit. 
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Table 4. DERIVED PARAMETERS OF THE AMMONIA STRUCTURES 



Component 


Ui,i) 


T^(2,2) 




T 

rot 


N(NH3) 


n(H2)^ 


M(H2) 










(pc) 


(K) 


(cm~^) 


(cm~^) 


(Mo) 


(Mo) 


C 


0.77 


0.49 


0.073 


24 


2.5 X 10^^ 


1 X 10^ 


16 


40 


NE 


0.39 


0.20 


0.029 


21 


1.0 X 10^5 


1 X 10^ 


1 


13 


S 


2.63 


<0.16 


0.050 


<10 


4.9 X 10^5 


3 X 10^ 


15 


13 



^Diameter estimated assuming a distance of 1.7 kpc. 
''Derived assuming [H2/NH3] = 10^ 



- 16- 

REFERENCES 

Aspin, C, & Geballe, T. R. 1992, A&A, 266, 219 
Beckwith, S. V. W., & Sargent, A. I. 1991, ApJ, 381, 250 
Beckwith, S. V. W., & Sargent, A. I. 1996, Nature, 383, 139 
Bonnell, I. A., & Bate, M. 2002, MNRAS, 336, 659 

Cesaroni, R., Felli, M., Testi, L., Walmsley, C. M., & Olmi, L. 1997, A&A, 325, 725 

Chernin, L. M., Masson, C. R., & Fuller, G. A. 1994, ApJ, 436, 741 

Devine, J. D., Bally, J., Reipurth, B., Shepherd, D., & Watson, A. 1999, AJ, 117, 2919 

Garay, G., & Lizano, S. 1999, PASP, 111, 1049 

Garay, G., Brooks, K., Mardones, D., & Norris, R. P. 2003, ApJ, 587, 739 

Girart, J. M., ct al. 1994, ApJ, 435, L145. 

Girart, J. M., Estalella, R., Viti, S., Williams, D. A., & Ho, P. T. P. 2001, ApJ, 562, L91 

Girart, J. M., Viti, S., Williams, D. A., Estalella, R., & Ho, P. T. P. 2002, A&A, 388, 1004 

Gomez, Y., Rodriguez, L. F., & Marti, J. 1995, ApJ, 453, 268 

Gyulbudaghian, A. L., Glushkov, Y. L., & Denisyuk, E. K. 1978, ApJ, 224, L137 

Heathcote, S., Reipurth, B., & Raga, A. C. 1998, AJ, 116, 1940 

Helmich, F. P., & van Dishoeck, E. F. 1997, A&AS, 124, 205 

Kurtz, S. 2000, RMxAAC, 9, 169 

Marti, J., Rodriguez, L. F., & Reipurth, B. 1993, ApJ, 416, 208 
Marti, J., Rodriguez, L. F., & Reipurth, B. 1995, ApJ, 449, 184 
Marti, J., Rodriguez, L. F., & Reipurth, B. 1998, ApJ, 502, 337 
Marti, J., Rodriguez, L. F., & Torrelles, J. M. 1999, A&A, 345, L5 
Reipurth, B., & Graham, J. A. 1988, A&A, 202, 219 



-17- 



Rodriguez, L. F., Moran, J. M., Dickinson, D. F., & Gyulbudaghian, A. L. 1978, ApJ, 226, 
115 

Rodriguez, L. F., Moran, J. M., Ho, P. T. P., & Gottlieb, W. 1980, ApJ, 235, 845 

Rodriguez, L. F., Marti, J., Canto, J., Moran, J. M., & Curiel, S. 1993, RexMexA&A, 25, 
23 

Rodriguez, L. F., Garay, G., Curiel, S., Ramirez, S., Torrelles, J. M., Gomez, Y., & Velazquez, 
A. 1994, ApJ, 430, L65 

Rodriguez, L. F., et al. 1998, Nature, 395, 355 

Rodriguez, L. F., Anglada, G., & Curiel, S. 1999, ApJS, 125 427 

Scoville, N. Z., & Kwan, J. 1976, ApJ, 206, 718 

Shepherd, D. S., Watson, A. M., Sargent, A. I., & Churchwell, E. 1998, ApJ, 507, 861 

Shepherd, D. S., Claussen, M. J., & Kurtz, S. E. 2001, Science, 292, 151 

Torrelles, J. M., Ho, P. T. P., Moran, J. M., Rodriguez, L. F., & Canto, J. 1986, ApJ, 307, 
787 

Torrelles, J. M., Gomez, J. F., Rodriguez, L. F., Curiel, S., Ho, P. T. P., & Garay, G. 1996, 
ApJ, 457, L107 

Trinidad, M. A., et al. 2003, ApJ, 589, 386 

Val'tss, I. E., Elhngsen, S. P., Slysh, V. I., Kalenskii, S. V., Otrupcek, R., & Larionov, G. 
M. 2000, MNRAS, 317, 315 

Viti, S., & Wilhams, D. A. 1999, MNRAS, 310, 517 

Viti, S., Girart, J. M., Garrod, R., Wilhams, D. A., & Estalella, R. 2003, A&A, 399, 187 

Wilner, D. J., Ho, P. T. P., & Rodriguez, L. F. 1996, ApJ470, L117 

Wilner, D. J., Reid, M. J., & Menten, K. M. 1999, ApJ, 513, 775 

Wilner, D. J., Ho, P. T. P., Kastner, J. H., & Rodriguez, L. F. 2000, ApJ, 534, LlOl 

Wiseman, J. J., & Ho, P.T.P. 1998, ApJ, 502, 676 

Yamashita, T., Suzuki, H., Kaifu, N., Tamura, M., Mountain, C. M., & Moore, T. J. T. 
1989, ApJ, 347, 894 



Zhang, Q., Hunter, T. R., & Sridharan, T. K. 1998, ApJ, 505, L151 

Zhang, Q., Hunter, T. R., Sridharan, T. K., & Ho, R T. R 2002, ApJ, 566, 982 



This preprint was prepared with the AAS IM^jX macros v5.0. 



-19- 




Fig. 1. — VLA image made at 7 mm of the exciting source of the HH 80-81 outflow. Contours 
are —3, 3, 4, 5, 6, and 8 x 1.3 mJy beam~^, the rms noise of the image. 
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Fig. 2. — Images of the BIMA continuum emission at 1.4 mm (top panel) and 3.5 mm 
(bottom panel). The contour levels are —3, 3, 4, 5, 6 8, 11, 14 and 17 x 0.038 and 
0.021 Jy beam~^, the rms noise of the 1.4 mm and 3.5 mm maps, respectively. The cross 
indicates the position of the thermal jet and the triangle marks the position of the H2O 
maser. The synthesized beams are shown in the bottom right corner of each panel. 
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Fig. 3. — Continuum cm and mm spectrum of the region associated with the exciting source. 
The dashed hne to the left is a linear fit to the four lower frequency data points, while the 
dashed line to the right is a linear fit to the three higher frequencies data points. 
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Fig. 4. — Image of the NH3 (1,1) velocity-integrated (from 10.0 to 14.3 km s~^) line emission 
from the surroundings of the exciting source of the HH 80-81 region. Contour levels are -3, 
3, 4, 5, 7, 9, 11 and 14 x 1.6 mJy beam~^, the rms noise of the image. This image is not 
corrected for the response of the primary beam and the S component appears weaker than 
it really is. The parameters given in Table 3 are corrected for the primary beam response. 
The ellipse indicates the position and alignment of the thermal jet, the triangle marks the 
H2O maser position taken from Gomez et al. (1995) and the small square to the south 
indicates the position of the radio continuum source VLA 13 corresponding to a knot of the 
jet (Marti et al. 1993). 
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Fig. 5. — Channel maps of the (1,1) (top) and (2,2) (bottom) inversion transition hnes of 
ammonia toward the HH 80-81 region. Contour levels are 4, 6, 8, 10, 15, 20, 30 and 50 x 
1.5 and 1.4 mJy bearn"^, respectively. The angular resolution is 5''1 x 3'.'0 at P.A.= —2°. 
The small cross indicates the position of the thermal jet and the triangle marks the position 
of the H2O maser. 
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Fig. 6. — Velocity integrated images of the new ammonia components; C, NE and S com- 
ponents. (Top and middle) C and NE clumps in the (1,1) and (2,2) ammonia emission 
integrated in a velocity range from 10.0 to 14.3 km s~^. The contour levels are -3, 3, 4, 
5, 7, 9, 11 and 14 x 1.6 and 1.5 mJy beam~^, the rms noise of the (1,1) and (2,2) maps 
respectively. The small cross indicates the position of the thermal jet, the triangle marks 
the position of the H2O masers and the large cross indicates the position and error of the 
methanol masers taken from Val'tts et al. (2000). (Bottom) Ammonia (1,1) emission of the 
S component integrated in a velocity range from 11.9 to 13.7 km s~^. The contour levels 
are -4, 4, 5, 6, 7, and 8 x 3.0 mJy beam~^, the rms noise of the map. This image was 
corrected for the response of the primary beam. The S component was detected only in the 
(1,1) ammonia transition. 
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Fig. 7. — Integrated spectra in the (1,1) and (2,2) lines of NH3 from the different molecular 
components in the HH 80-81 region. The small vertical bars indicate the expected velocities 
of the main and satellite (1,1) ammonia lines. 
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Fig. 8. — Bottom panels. Images of the BIMA integrated emission of the SO 65-44 (left 
panel), HCO+ 1-0 (central panel), HON 1-0 (right panel) lines. The SO map shows the 
integrated emission in a velocity range from 10.0 to 14.3 km s~^. The HCO"'" and HON 
maps show the integrated emission in a velocity range from 9.5 to 15.7 km s~^. The contour 
levels are —2, 2, 3, 4, 5, 6, 8, 10, 12 and 14 x 3.2, 0.15, 0.13 Jy beam^-*^, the rms noise of 
the SO, HCO+ and HCN maps, respectively. The small cross indicates the position of the 
exciting source of the thermal jet and the triangle marks the position of the H2O maser. 
The large cross indicates the position and error of the methanol masers taken from Val'tts 
et al. (2000). The synthesized beams are shown in the bottom right corner of each panel. 
Top panels. Spectra of the SO 65-44 (left panel), HCO+ 1-0 (central panel), HCN 1-0 (right 
panel) lines. The SO spectrum was taken at the position of the thermal radio jet and divided 
by a factor of 5 in order to put it at the same scale than the HCO"'' and HCN spectra that are 
weaker than the SO. For the HCO"*" and HCN the spectra was obtained from the averaged 
emission within the dashed box shown in the HCN map. 
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Fig. 9. — Position-velocity diagrams of the NH3 (1,1) emission from the C component along 
the major axis direction at P.A.=14.7° (top) and along the minor axis direction P.A. =104.7° 
passing through a(2000)= 18^ 19™ 12f5; 5(2000)= -20° 47' 24'.'9. The contour levels are 
-4, 4, 6, 8, 10, 15, 20, 30, and 40 x 1.4 mJy beam-i. 
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Fig. 10. — Ammonia spectra in the (1,1) and (2,2) lines at the peak position of the HH 80-81 
thermal jet: q;(2000)= 18'* 19"* 12f509; 5(2000)= -20° 47' 27'/42. Only weak and broad (2,2) 
emission was detected. 



